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"Glauber" models are used to calculate geometric quantities in the initial state of heavy ion collisions, 
such as impact parameter, number of participating nucleons and initial eccentricity. The four RHIC 
experiments have different methods for Glauber Model calculations, leading to similar results for 
various geometric observables. In this document, we describe an implementation of the Monte Carlo 
based Glauber Model calculation used by the PIfOBOS experiment. The assumptions that go in 
the calculation are described. A user's guide is provided for running various calculations. 



I. INTRODUCTION 



A. Makeup of Nuclei 



In heavy-ion collisions, initial geometric quantities 
such as impact parameter and shape of the collision re- 
gion cannot be directly determined experimentally. How- 
ever, it is possible to relate the number of observed par- 
ticles and number of spectator neutrons to the centrality 
of the collision. Using the percentile centrality of a colli- 
sion, the initial geometric configuration can be estimated 
with models of the contents of a typical nucleus. 

These models fall in two main classes. (For a recent 
review, see Ref. [T].) In the so called "optical" Glauber 
calculations, a smooth matter density is assumed, typi- 
cally described by a Fermi distribution in the radial di- 
rection and uniform over solid angle. In the Monte Carlo 
based models, individual nucleons are stochastically dis- 
tributed event-by-event and collision properties are cal- 
culated by averaging over multiple events. As discussed 
in Ref. [1] and Ref. [5] , these two type of models lead to 
mostly similar results for simple quantities such as the 
number of participating nucleons (A^part) and impact pa- 
rameters (6), but give different results in quantities where 
event-by-event fluctuations are significant, such as par- 
ticipant frame eccentricity (epart)- 

In this paper, we discuss in detail the Monte Carlo 
Glauber calculation implemented by PHOBOS. In sec- 
tion |ll] the method is outlined and the assumptions that 
go into the calculation arc introduced. In section |III| 
we discuss the implementation and the tutorial functions 
provided. 



II. THE MODEL 



The Monte Carlo Glauber Model calculation is per- 
formed in two steps. At first, the nucleon positions in 
each nucleus are stochastically determined. Then, the 
two nuclei are "collided", assuming the nucleons travel 
in a straight line along the beam axis (eikonal approxi- 
mation) such that nucleons are tagged as wounded (par- 
ticipating) or spectator. 



The position of each nucleon in the nucleus is deter- 
mined according to a probability density function. In 
a quantum mechanical picture, the probability density 
function can be thought of as the single-particle proba- 
bility density and the position as the result of a position 
measurement. In the determination of the nucleon posi- 
tions in a given nucleus, it is possible to require a mini- 
mum inter-nucleon separation (dmin) between the centers 
of the nucleons. 

The probability distribution is typically taken to be 
uniform in azimuthal and polar angles. The radial prob- 
ability function is modeled from nuclear charge densi- 
ties extracted in low-energy electron scattering experi- 
ments [5]. The nuclear charge density is usually param- 
eterized by a Fermi distribution with three parameters: 



l + w{r/R)^ 



(1) 



where po is the nucleon density, R is the nuclear radius, 
a is the skin depth and w corresponds to deviations from 
a spherical shape. The overall normalization (po) is not 
relevant for this calculation. Values of the other param- 
eters used for different nuclei are listed in Table HI 

Two exceptions are the deuteron (^H) and sulfur ("^^S) 
nuclei. For sulfur, a three parameter Gaussian form is 
used: 



Pir) ^ Po 



1 + wir/Rf 
1 + exp{ - 



(2) 



The values of R, a and w for sulfur are also given in 
Table |Tj For deuteron, three options are supported: 

1. The three parameter Fermi distribution can be 
used, with the values given in Table |l] 

2. The Hulthen form can be used: 

p[r) = Po , (3) 



where a = 0.457 fm^^ and b ^ 2.35 fm^^ g] |5]. 

The proton can be randomly placed using the 
Hulthen form given above and the neutron can be 
placed opposite to it. 
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Nucleus 


R [fm] 


a [fm] 


w [fm] 




0.01 


0.5882 





16q 


2.608 


0.513 


-0.51 




3.34 


0.580 


-0.233 


32g 


2.54 


2.191 


0.16 


«Ca 


3.766 


0.586 


-0.161 




4.309 


0.517 


-0.1308 




4.2 


0.596 





1S6W 


6.58 


0.480 







6.38 


0.535 





207p|^a 


6.62 


0.546 





238 u 


6.81 


0.6 






"These values are also used for ^"^Pb for which Fermi parameters 
are not available. It has been noted that Bcsscl-Fouricr coefficients 
for the two nuclei are similar [S]- 

TABLE L Nuclear charge density parameters for different nu- 
clei, taken from Ref. [3]. 



It should be noted that the 3'''* option was used in PHO- 
BOS analyses. 
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B. Collision Process 

The impact parameter of the collision is chosen ran- 
domly from a distribution dN/db cx 6 up to some large 
maximum 6max with femax — 20 fm> 2Ra ■ The centers of 
the nuclei are calculated and shifted to (—5/2,0,0) and 
(6/2,0,0) ^. It is assumed that the nucleons move along 
a straight trajectory along the beam axis. (The longitu- 
dinal coordinate does not play a role in the calculation.) 

The inelastic nuclcon-nucleon cross section (ctn^), 
which is only a function of the collision energy is ex- 
tracted from p-l-p collisions. At the top RHIC energy of 
^Snn = 200 GeV, ctnn — 42 mb, while at the LHC it is 
expected to be around CTnn = 72 mb (with large uncer- 
tainty from the unknown elastic cross section) . The "ball 
diameter" is defined as: 

D = yJoNNh- (4) 

Two nucleons from different nuclei are assumed to collide 
if their relative transverse distance is less than the ball 
diameter. If no such nucleon-nucleon collision is regis- 
tered for any pair of nucleons, then no nucleus-nucleus 
collision occurred. Counters for determination of the to- 
tal (geometric) cross section are updated accordingly. 
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FIG. 1: Typical events for Cu-|-Cu (top panel), Au-I-Au (mid- 
dle panel), and Pb+Pb (lower panel) collisions, the first two 
performed at RHIC energies and the latter at the LHC. 
Wounded nucleons (participants) are indicated as solid cir- 
cles, while spectators are dotted circles. 



III. USERS' GUIDE 



Throughout the paper, the reaction plane, defined by the impact 
parameter and the beam direction, is given by the x- and z-axes, 
while the transverse plane is given by the x- and i/-axes. 



The PHOBOS Glauber MC code works 
within the ROOT framework (ROOT 4.00/08 or 
higher [B]). The code is contained in the macro 
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runglauber_vX. Y.C [7] (Latest version is 1.0.). Three 
classes, TGlauNucleon, TGlauNucleus and TGlauberMC 
and two example functions runAndSaveNtupleO and 
runAndSaveNucleons arc defined in the provided 
macro. While the functionality is essentially complete 
for known applications of the Glauber approach, users 
are encouraged to write their own functions to access 
results of the Glauber simulation or to modify the code: 

• TGlauNucleon is used to store information about a 
single nucleon. The stored quantities are the posi- 
tion of the nucleon, the number of binary collisions 
that the nucleon has had and which nucleus the nu- 
cleon is in, "A" or "B". For every simulated event, 
the user can obtain an array containing all nucle- 
ons (via TGlauberMC: :GetNucleons()). 

• TGlauNucleus is used to generate and store infor- 
mation about a single nucleus. The user is not ex- 
pected to interact with this class. 

• TGlauberMC is the main steering class used to gen- 
erate events and calculate event-by-event quantities 
such as the number of participating nucleons. 

The steering class TGlauberMC has one constructor 

TGlauberMC: : TGlauberMC (Text_t* NA, 

Text_t* NB, 
Double_t xsect) 

where NA and NB are the names of the colliding nuclei and 
xsect is the nucleon-nucleon cross section given in mb. 
The defined nuclei names are: "p" , "d" , "dhh" , "dh" , 
"O", "Si", "S", "Ca", "Cu", "W", "Pb", "Au", "Ni" and 
"U" (see Table For deuteron, the names "d", "dhh" 
and "dh" correspond to the three options described in 



section II A respectively. Units are generally given in fm 
for distances, while in mb for cross sections. 



A. Running the Code 

To generate Au-fAu collisions at ^s„„ = 
200 GeV (^Snn — 42 mb) one would construct a 
TGlauberMC object by issuing the commands: 

root [0] .L runglauber_X. Y.C+ 

root [1] TGlauberMC glauber("Au","Au", 42) ; 

where the first ROOT command compiles, links and loads 
the compiled macro ^ including the Glauber code as ex- 
plained in chapter 2 of the ROOT users' guide. 

Events can be generated interactively using the two 
functions 




part 




FIG. 2: Distributions of A^part and A^coii for 10k events for 
Cu-I-Cu and Au-|-Au at RHIC, and Pb-I-Pb at the LHC. 



• TGlauberMC: :NextEvent(Double_t bgen) , which 
is used to run an event at a specified impact pa- 
rameter, or over a range of impact parameters (if 
bge n=-l, the default value) as described in sec- 
tion 



IIBl 



• TGlauberMC: :Run(lnt_t nevents) which is used 
to run a large event sample by invoking NextEvent 
many times. 

Other important public member functions are: 

• TGlauberMC: :SetMinDistance(Double_t d), 

which is used to set minimum nucleon seperation 
within a nucleus, dmin (default is 0.4 fm) 

• TGlauberMC: :SetBmin(Double_t bmin) and 
TGlauberMC: :SetBmax(Double_t bmax), which 
can be used to set the range of impact parameter 
values generated in Run(). 

• TGlauberMC: :GetTotXSect() which returns the 
total nucleus-nucleus cross section, calculated when 
the function RunO is called. 



^ Note that you must replace X.Y with the current version number 
of the code, for example 1.1. 



• TGlauberMC: :Draw() which 
event in the current pad. 



draws the current 
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B. Example functions 

Two example functions are provided to demonstrate 
how to run the model. 

runAndSaveNtupleO generates a number of Monte 
Carlo events and saves some event-by-event quantities. 
It takes as parameters, the number of events to be gener- 
ated, the collision system, the nucleon-nucleon cross sec- 
tion, the minimum separation distance and the output 
file name. It creates and stores an ntuple in the output 
file with the following event-by-event quantities: 

• Npart: Number of participating nucleons. 

• Ncoll: Number of binary collisions. 

• B: Generated impact parameter. 

• MeanX: Mean of x for wounded nucleons, {x}. 

• MeanY: Mean of y for wounded nucleons, {y}. 

• MeajiX2: Mean of for wounded nucleons, {x^}. 

• MeanY2: Mean of y^ for wounded nucleons, {y^}- 

• MeanXY: Mean of xy for wounded nucleons, {xy}. 

• VarX: Variance of x for wounded nucleons, cr^. 



• VarY: Variance of y for wounded nucleons, cr^. 



• VarXY: Covariance of x and y for wounded nucleons, 
a.^y = {xy} - {x}{y}. 

• MeanXSystem: Mean of x for all nucleons. 

• MeanYSystem: Mean of y for all nucleons. 

• MeanXA: Mean of x for nucleons in nucleus A. 

• MeanYA: Mean of y for nucleons in nucleus A. 

• MeanXB: Mean of x for nucleons in nucleus B. 
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FIG. 3: eRp (open symbols) and tpart (closed symbols) as a 
function of A^part for Cu+Cu and Au+Au collisions at RHIC 
and Pb-|-Pb collisions at the LHC. 



C. Sample Results 

As an example application of this code, 10k events 
were generated for Cu-|-Cu and Au+Au at RHIC ener- 
gies (^Snn = 42 mb), and Pb+Pb at LHC beam en- 
ergy i\/sm ~ 72 mb) using the runAndSaveNtupleO 
function. The resulting ntuples were used to plot the 
distributions of A^part and A'coii, shown in Fig. |2] Using 
the event-by-event quantities, one can construct combi- 
nations of moments like [2] : 

• Reaction-plane eccentricity crp 

VarY - VarX 



VarY + VarX 
• Participant eccentricity epart 



(5) 



'(VarY - VarX) 2 + 4VarXY2 
VarY + VarX 

which are shown in Fig. [3] for the different systems. 



(6) 



• MeanYB: Mean of y for nucleons in nucleus B. 

It is important to note that for each of these event- 
by-event quantities a "getter" function is implemented 
providing the users the option to write their own event 
loop (using TGlauberMC: :NextEvent () . 

The function runAndSaveNucleons () generates a 
number of Monte Carlo events and saves an array of 
TGlauNucleon objects for each event. It is also possi- 
ble to use this function to print out the values stored 
in the nucleons by setting the verbosity parameter. The 
function takes as parameters the number of events to be 
generated, the collision system, the nucleon-nucleon cross 
section, the minimum separation distance, the verbosity 
flag and the output flle name. 



IV. CONCLUSION 

This work has described the PHOBOS implementation 
of the "Glauber Model" commonly used by heavy ion 
physics experiments to study the initial state configura- 
tions of nuclear matter. The code, accessible online, can 
be used within user code or in a standalone mode allow- 
ing analysis of various distributions (e.g. Npart, Ncoii, 
b, £part)- The authors welcome comments on the code 
and suggestions on how to make it more useful to both 
experimentalists and theorists. 
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